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DERIVATION  OF  EQUATIONS  FOR  CONVERTING  FROM  GEODETIC 
COORDINATES  TO  GEOCENTRIC  COORDINATES 

F.  T.  Heuring 

In  the  A.P..L.  " orbit  computation  programs,  the  TRANET  Tracking 
Stations  are  specified  in  a  geocentric  coordinate  system,  whereas, 
particular  positions  (such  as  a  TRANET  Tracking  Site)  over  the  Earth  are 
expressed  initially  in  a  geodetic  coordinate  system.  In  order  to  acquire 
geocentric  coordinates  from  a  given  set  of  geodetic  coordinates  a  set  of 
transformation  equations  were  derived. 

Section  I  will  define  the  notation,  and  Section  II  will  embody 
the  derivation  of  the  transformation  equations. 

* 

I.  Notation 

Let ;  , 

9^^  =  geodetic  latitude  of  i-th  tracking  site  in  its  local  datum, 

=  geodetic  longitude  of  i-th  tracking  site  in  its  local  datum, 

h^  =  elevation  of  i-th  tracking  site  above  (below)  geoid, 

=  geoidal  height  of  i-th  tracking  site  in  its  local  datum, 

=  deflection  in  meridian  at  i-th  tracking  site, 

T|^  =  deflection  in  prime  vertical  at  i-th  tracking  site, 

a^  =  equatorial  radius  of  the  local  datum  spheroid  of  the  i-th 
tracking  site,  scaled  by  R^, 

b^  =  polar  radius  of  the  local  datum  spheroid  on  the  i-th  tracking 

site,  scaled  by  R  , 

^  •  o 

See  References  1  and  2  for  definition  of  geodetic,  datum,  etc. 
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X 

Ax^,  Ay^,  Az^ 


cartesian  coordinates,  scaled  ty  Rq,  on  i-th  datum  spheroid 
as  specified  by  tracking  site  cp^  and  (cartesian  origin 
identical  to  i-th  datum  origin), 

cartesian  coordinates,  scaled  by  R^,  on  geoid  as  specified  by 
tracking  site  H,  (cartesian  origin  identical  to  i-th  datum 
origin), 

cartesian  coordinates,  scaled  by  Rq,  of  tracking  site  on 
earth’s  surface,  (cartesian  origin  identical  to  i-th' datum 
origin). 


center  of  spheroid  of  the  i-th  tracking  site  datirni  in  the 
A.P.L.  Datum,  scaled  by  R^, 


equatorial  radius  of  A.P.L.  Datum  spheroid, 


cartesian  coordinates  of  tracking  site  in  A.P.L.  geocentric 
coordinates,  scaled  by  R^, 

radius  of  i-th  tracking  site  in  A.P.L.  geocentric  coordinates, 

scaled  by  R  , 
o 

latitude  of  i-th  tracking  site  in  A.P.L.  geocentric  coordinates. 


longitude  of  i-th  tracking  site  in  A.P.L.  geocentric  coordinates. 


(x 


2 

y.  )  • 


II.  Derivation 

A.  Given  ^  ^  conversion  to  ,  y^  ,  and  is  as 

follows.  Using  the  equation  for  an  ellipse 
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in  particular  the  ellipse  is  a  meridional  plane  of  the  i-th  datum; 
differentiate  z„  with  respect  to 

Ct.  (jt  . 


But  (see  Figure  lA), 


tan  cpg 


from  which  by  algebraic  manipulation  .(^’igure  IB), 


(1  +  (-^f  tan^  9p 


after  which, 


cos 


sin 


tan  cp^j  . 


A 


y 


\ 

f 


(1) 


pigur*  1A 

1 

♦o"  i>Q. 


Meridian  Plan*  in  i-th  Datum 
^  *G: 


=  cat  ^Gj 


=*  tan  e  — 


^G: 
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Figure  2 

Diggrgm  of  tho  deflections  of  the  vertical  (^j  and  r/j )  and  the  associated 
quantities  necessary  to  acquire  the  cartesian  coordinates  on  the  geoid 
from  earth  surface  cartesian  coordinates. 
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D.  Let  us  simplify  by  expanding  small  quantities.  Assume: 


\  i  30"  (of  arc); 


1°  <  |cp(j  1  <  89°; 

i 


and  only  take  quantities  of  magnitude  T]^  and  to  second  order. 


cos  =  1  -  — — ,  sin 


cos  Tlj^  i  1  -  ,  sin  Tlj^  i 


cos  =  1  -  — 5- 


thus, 


c 

(9q  +  1^)  =  Ceos  cpg  (1  -  sin  cp^  T 


2  2  2  2 
=  (1  -  — )  cos  cp  +  §  sin  cp„  ■ 

i  i 


-  2  §^(l  >;  — )  sin  cos  9^ 

1  i 


2  2 

=  cos  q)p  -  sin  2(p„  -  cos  29-  +  3rd  order 


From  a  personal  communication  with  Mr.  L.  Simmons,  U.S.C.  and  G.S., 
deflection  of  30"  exist  but  are  in  general  uncommon. 
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F.  The  cartesian  coordinates  in  the  A.P.L.  geocentric  system  are: 


X  =  x_,  +  Ax. 

Cl  ^1 


^c,  =  ^E.  +  *^1 
i  1 


^c.  = 

1  1 


where  Ax^,  ^^i^  ^^i  second  order,  at  "best. 


H.  The  cylindrical  coordinates  (z  ,  Q  ,  X  )  in  the  A.P.L. 

'^i  '^i  '^i 

geocentric  system  are: 


i 

~  ~2 —  ^^i  ^i^  '’’g  ^i  ^i  '’’g  '^■^i  *  (l3) 


-2  2  2 
C„  =  Xc  +  yc 
°i  i  i 


After  some  algebraic  manipulation  and  using  the  binominal  expansion 


=  Cg  ^^i  ^1^  'Pg  ^^i  ^G  ^^i  ^G  "  ^1  ^i  '^’g 


+  +  hj^)  cos  '  —pi —  (Ax^  cos  +  Ay^  sin  \q  ) 

i  u .  1  i 


.+  3rd  order. 
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Figure  3A  Diagram  shewing  means  of  determining  when  hj  =  0. 


APL  DATUM 
ORIGIN 


Figure  3B 


Diagram  showing  datarmination  of  whan  hj  4  0- 
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and  it  follows  slmilafly 


From  equation  (15)  and  Figure  3B, 


i  i  1 

=  Xp  +  — p - p -  [Ax  sin  (\»  +  p)  -  Ay  cos  (\  +,  p)] 

^i  ^G.  ^G^  ^  ^i  ^  ‘^i 

2 

Assuming  cos  p  =  1 - 1 ,  sin  p  =  p, 


X^  =  [h^  Tl^  -  (Ax^  sin  -  Ay^  cos  X^  )]  +  3rd  order. 


Equations  (13),  (l^)^  and  (16)  are  the  cylindrical  coordinates  z 

^i 

C  ,  X.  in  the  A.P.L.  Earth  fixed  coordinate  system  expressed  as  a  function 
i  i 

of  the  geodetic  coordinates  of  a  tracking  station. 


(16) 
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